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We study theoretically the equilibrium phase behaviour of a mixture of polydisperse hard-sphere 
colloids and monodisperse polymers, modelled using the Asakura-Oosawa model within the free 
volume approximation of Lekkerkerker et al. We compute full phase diagrams in the plane of colloid 
and polymer volume fractions, using the moment free energy method. The intricate features of phase 
separation in pure polydisperse colloids combine with the appearance of polymer-induced gas-liquid 
coexistence to give a rich variety of phase diagram topologies as the polymer-colloid size ratio £ 
and the colloid polydispersity 5 are varied. Quantitatively, we find that polydispersity disfavours 
fluid-solid against gas-liquid separation, causing a substantial lowering of the threshold value £ c 
above which stable two-phase gas-liquid coexistence appears. Phase splits involving two or more 
solids can occur already at low colloid concentration, where they may be kinetically accessible. 
We also analyse the strength of colloidal size fractionation. When a solid phase separates from a 
fluid, its polydispersity is reduced most strongly if the phase separation takes place at low colloid 
concentration and high polymer concentration, in agreement with experimental observations. For 
fractionation in gas- liquid coexistence we likewise find good agreement with experiment, as well as 
with perturbative theories for near-monodisperse systems. 
PACS numbers: 64.70Fx, 68.35.Rh 



I. INTRODUCTION 

Since the 1980s phase separation phenomena in col- 
loidal suspensions that are induced by the addition of 
polymers have gradually gained interest both from a tech- 
nological and scientific point of view. The equilibrium 
phase behaviour of a pure colloidal hard sphere system 
is dictated exclusively by entropic effects, and these can 
cause separation into fluid and solid phases. The term 
fluid is used to emphasise that there are no separate gas 
and liquid phases in such a system. However, such phases 
can appear with the addition of polymer, because of 
the effective colloid-colloid attraction which this induces. 
Our aim in this paper is to study the effect in which 
colloid polydispersity, i.e. a continuous spread of colloid 
sizes, further affects the phase behaviour. Because poly- 
dispersity causes complex phase equilibria even for pure 
hard spheres [1, 2], with for example the coexistence of 
a fluid with several solids, one expects very rich phase 
behaviour when polymer is added. 

A theoretical model of phase separation in a mixture 
of colloidal and polymer particles was first advanced by 
Asakura and Oosawa (AO) [3] and extended by Vrij [4]. 
In order to describe the interactions in a solution with 
non-adsorbing polymers, they proposed to model each 
polymer chain as a sphere with diameter er p equal to 
twice the radius of gyration. These "polymer spheres" 
are then assumed to be able to interpenetrate freely with 
each other, which is a reasonable assumption for poly- 
mers near their 0-point where the chains obey random 
walk statistics. On the other hand, the polymer chains 
are completely excluded from the space occupied by the 
colloidal particles, corresponding to a hard interaction 
between polymers and colloids. 

Overall, the AO model is effectively a mixture of hard 
spheres, but with extreme non-additivity in the polymer- 



polymer interaction which has range zero rather than 
(T p . Explicitly, the AO interparticle potentials can be 



expressed as follows: 



v cc (r) = 
v cp (r) = 



if r < i(cr c +<t' c ) 



00 

otherwise 



if r < §(cr c + o-p) 



oo 

otherwise 



pp 



(r) = 



Here r is the distance between the centres of mass of 
the particles considered. The subscript "cc" refers to an 
interaction between two colloidal particles; in a polydis- 
perse system these can have different diameters er c and 
o~' c . Similarly "cp" and "pp" denote colloid-polymer and 
polymer-polymer interactions, respectively. 

Even though all the interactions are still hard within 
this model, the entropy of the polymers induces an ef- 
fective attraction between the colloids. The centres of 
mass of the polymers are excluded from a spherical "ex- 
clusion zone" of width er p /2 around the colloids. When 
the exclusion zones of two or more colloids overlap, the 
polymers cannot access the region between the colloids. 
The osmotic pressure of the polymers is then unbalanced 
and creates an attractive depletion interaction between 
the colloids. Its range is set by the polymer diameter tr p , 
and its strength at contact is found to be proportional 
to (TpCr c times the osmotic pressure of the polymers. It is 
this depletion interaction which can cause the suspension 
to separate into colloid-poor and colloid-rich phases, re- 
sulting in phase equilibria with solid, liquid or gas phases 
depending on conditions. These striking effects have been 
extensively investigated experimentally [5-9] and the AO 
model has been shown to give a good description for 
suspensions of sterically stabilised colloidal particles im- 
mersed in a solvent with polymers [10, 11]. However, ex- 
perimental systems are always polydisperse, containing a 
spread of colloid sizes as well as polymer chain lengths. 
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FIG. 1: Phase diagram sketch for colloid-polymer mixtures 
with three different values of the size ratio £. For £ < £ c 
the depletion interaction is weak and the only effect is the 
widening of the (F)luid-(S)olid coexistence region; at £ > £ c 
(right) the longer polymers create a sufficiently long-range 
interaction that induces the formation of (G)as and (L)iquid 
phases and hence of a critical point, marked with a circle. The 
crossover value £ ~ £ c (middle) is the value where a region of 
gas-liquid coexistence first appears. 



As reviewed below, our theoretical understanding of the 
effects of this polydispersity on the phase behaviour re- 
mains very limited, fdeally, one would therefore like to be 
able to study theoretically the phase behaviour of a mix- 
ture of colloids and polymers that are both polydisperse. 
This is a very challenging problem. As a first step we 
therefore focus in this paper on the case where only the 
colloids are treated as polydisperse. This already turns 
out to result in very rich phase behaviour. 

We first review some of the previous work carried out 
on colloid-polymer mixtures. Most theoretical studies 
have focused on the case where the colloids are monodis- 
perse and the polymers are either monodisperse or poly- 
disperse. Theoretical work by Gast et al [12], Lekkerk- 
erker et al [10] and Dijkstra et al [11] for a monodisperse 
mixture predicts that the topology of the phase diagram 
depends crucially on the range a p of the depletion in- 
teraction, or more precisely the ratio £ = cr p /<T c of the 
polymer and colloid diameters. This is sketched in Fig. 1, 
which shows phase diagrams in terms of colloid volume 
fraction <f> c and effective polymer volume fraction 4> p . The 
latter is defined as the polymer concentration p p times 
the volume (7r/6)(7p of each polymer and can be greater 
than 1 because the polymers can overlap. 

To understand these phase diagrams, one can imag- 
ine gradually adding polymers to a pure (monodisperse) 
colloidal hard-sphere system which exhibits coexistence 
between a fluid and a crystalline solid with volume frac- 
tions of around 50% and 55%, respectively [13-15]. For 
size ratios £ below a certain crossover value £ c the addi- 
tion of polymer only has the effect of widening the fluid- 
solid coexistence region; see Fig. 1 (left). This widening 
is gradual at first but then becomes very pronounced. It 
has been argued that, when viewed in terms of the colloid 
system with its effective polymer-induced attraction, this 
can be understood as the crossover from a regime with 
dominant entropic effects to an "energetic fluid" regime 
where the strength at contact of the attractive interac- 
tion is the key parameter [16]. 



For sufficiently large polymers, on the other hand 
(£ > £c), the addition of polymeric particles leads to a 
new phase diagram topology where gas-liquid (and there- 
fore gas-liquid-solid) coexistence can occur, as shown in 
Fig. 1 (right). The threshold value £ c of the size ratio is 
the one where a gas-liquid region first appears; see Fig. 1 
(middle). The intuitive explanation for the different be- 
haviour at small and large £ is that if the interaction 
is generated by sufficiently long polymers, £ > £ c , then 
neighbouring colloids can still be within the range of the 
attraction even when they are in a disordered state, and 
so a thermodynamically stable liquid phase can form. 

The actual crossover value £ c varies slightly depend- 
ing on the theory considered but is found to lie in the 
range £ c w 0.3 . . . 0.6. E.g. in [11] a colloid system with 
the effective AO interaction was simulated with the re- 
sult 0.4 < £ c < 0.6; this result is approximate because in 
this range of £ the effective interaction potential for the 
colloids also contains many-body terms. Lekkerkerker et 
al [10], by the use of a van der Waals-type approximation 
(see below) identified £ c w 0.32, and Gast and collabora- 
tors [12], with the aid of perturbation theory, estimated 
£ c to be « 1/3. On the experimental side Ilett et al [5] 
found instead £ c w 0.25 for a suspension with polydisper- 
sity S of order 0.05. One of the goals of our study will be 
to identify the effects of colloidal polydispersity on the 
threshold value £ c . 

Initial attempts to understand the influence of poly- 
dispersity on colloidal-polymer mixtures have been made 
in recent years, but with a focus on size-polydispersity 
in the polymers [17-20] which causes quantitative rather 
than qualitative changes in the phase diagrams. To the 
best of our knowledge the only studies considering poly- 
disperse colloids investigate the liquid-gas phase sepa- 
ration in the limit of near-monodisperse colloids, where 
perturbative theories can be applied [9, 21, 22]. In this 
regime it turns out that the gas-liquid coexistence region 
widens as the colloidal spheres becomes more polydis- 
perse [22] . This statement applies when the osmotic pres- 
sure of the polymer - or equivalently its chemical poten- 
tial - is imposed, e.g. by connecting the colloid-polymer 
mixture to a large polymer reservoir. This corresponds 
to fixing the effective interaction between the colloids. 
More recent experimental work [9] has focused on the 
strength of size fractionation between coexisting phases 
and its scaling with the overall colloid polydispersity; we 
return to a comparison with this study in Sec. VB. 

We provide here a theoretical analysis of phase equi- 
libria in colloid-polymer mixtures with significant size- 
polydispersity in the colloids. The paper is structured 
as follows. In Sec. II we derive the free energy of the 
polydisperse AO model, within the free volume approx- 
imation employed by Lekkerkerker et al [10]. This free 
energy has a truncatable [23] structure, allowing us to 
solve the resulting phase equilibrium equations by using 
the moment free energy (MFE) method as explained in 
Sec. III. In Sec. IV we study the phase diagram topolo- 
gies that arise, while Sec. V is devoted to a quantitative 
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analysis of the effects of polydispersity and polymer size 
on the phase diagrams, and of colloidal size fractionation 
between coexisting phases. Conclusions and an outlook 
towards future work are given in Sec. VI. 



II. FREE VOLUME THEORY 

Even though in the AO model the polymer-polymer 
interactions are ideal, a direct calculation of the parti- 
tion function is a difficult task because one needs to keep 
track of the overlaps of all exclusion zones around the 
colloids; the remaining free volume accessible to the poly- 
mers thus depends in a complicated way on the configu- 
ration of the colloids. Even if only pairwise overlaps are 
considered one still effectively has a colloid system with 
a depletion interaction potential of finite range, which 
cannot be solved exactly. Progress can however be made 
by using a van der Waals or free volume approximation 
which replaces the free volume by its average over colloid 
configurations. This will lead us to an explicit expression 
for the Hclmholtz free energy of our system as a sum of 
two contributions, one corresponding to the pure hard- 
sphere system and one describing the average interaction 
between polymers and colloids. 

In what follows we will derive the free volume ap- 
proximation to the free energy for a system where both 
the polymers and colloids are polydisperse. Allowing for 
polydispersity in the polymers adds almost no extra work 
and facilitates comparison with other approaches. The 
case of monodisperse polymers which we later study nu- 
merically is easily recovered as a special case. 

We assume initially that we have a discrete mixture of 
colloid and polymer species and take the polydisperse 
limit at the end. The colloid species are labelled by 
s, with {-/Vcs} the number of particles of each species 
and a cs their diameter; similarly {N pt } and a p t give the 
number of particles and diameter for each of the poly- 
mer species. Because all interactions in the AO model 
are hard, temperature becomes an irrelevant scale factor 
which only sets the energy scale. We therefore measure 
all energies in units of T = 1/(3 and set ks to one. The 
spatial coordinates of the particles are denoted by Ri for 
the colloids and rj for the polymers. The canonical par- 
tition function is then 



Z = 



n 



A, 



pt I dR dr 



.! N pt \ 



x cxp 



(1) 



where v cc (i,i') stands for v cc (\Ri — Ri>\) and similarly 
for v cp (i,j). The integrations dR, dr extend over all col- 
loid and polymer coordinates, and A cs and A p t are the de 
Broglie wavelengths for each species. The natural ther- 
modynamic potential to work with, in this context, is the 
semi-grand canonical potential. For this the system can 



be thought of as connected to a polymer reservoir via 
a membrane impermeable to colloids but permeable to 
polymers, so that the polymer particle numbers {N pt } 
can range over all possible values. The resulting semi- 
grand partition function is given by 

s({AU,v,{ Mpt })= E l[e^ N »*z({N cs },{N pt },v) 

{N pt } t 

(2) 

where the chemical potentials fi p t of the polymers are 
fixed by the reservoir. In (2) we can separate the ideal 
contributions from the excess part: 

s = *£.E /^«p(-E««(m')) 

{N pt } J \ »<»' J 

«n£(/ 



dr 



Ejj «cp(*>j) 



(3) 



where Z® s is the canonical partition function of an ideal 
mixture of colloids, Z° s = V N * ]J S Ks 3Ncs and N c is 
the total number of colloids. The quantity e Mpt /X pt iden- 
tifies the thermodynamic activity of the polymer species 
t, which for our ideal polymers is just its density p pt 
in the reservoir. The integration over Vjm, which rep- 
resents the position of any one of the polymers from 
species t, gives the fraction of volume available to a hard 
sphere of diameter a pt if the colloids are in a configura- 
tion {Ri}, a({Ri},a pt ) = V" 1 j dr j{t) e-^^\ R ^i\\ 
We are thus left with a configurational integral over the 
hard spheres: 

xexpJvE^^W)) ' W 

Following Widom [24] this can be reinterpreted as an av- 
erage ( )h s over the Boltzmann distribution of the pure 
colloidal hard-sphere system, with excess partition func- 



tion Zg: 



— ' r^O ycx 



cxp (v^2 p pt Q ({- R i}' <7 pt)^ 



(5) 



hs 



We note that the second exponential in (4) is the effec- 
tive colloid-colloid interaction generated by the polymers. 
Beyond trivial zero- and one-body terms it always con- 
tains a two-body contribution which is the attractive de- 
pletion interaction described above. For sufficiently large 
polymers (£ > 0.154) also interactions between three or 
more colloids appear because more than two exclusion 
zones can overlap simultaneously. 

The average in (5) cannot be evaluated exactly. To 
make progress, we follow the van der Waals approach or 
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"free volume theory" of Lekkerkerker et al [10]. This con- 
sists of "moving" the average ( )h s inside the exponential 
to give 



S = i^£exp (vJ2pU W{Ri},a pt )) 



lis 



(6) 



This approximation is exact in the limit of small poly- 
mer activities p pt — > 0; otherwise it can be viewed as 
identifying a lower bound on S, or an upper bound on 
the thermodynamic potential £1 = — InS. The latter be- 
comes 



n = F hs -Vj2pU (a({Ri},a pt )) 



hs 



(7) 



where Fh s is the Hclmholtz free energy of the pure colloid 
system. The average free volume (a({Ri}, o~ p t)) hs which 
appears here is the probability of being able to insert a 
hard particle of diameter a p t into the pure colloid sys- 
tem. But from the Widom insertion principle [24] this is 
just exp[— /i^(c7 pt )] where P^(<t) is the excess chemical 
potential of a colloid particle of diameter a, in a system 
without polymer. The semi-grand potential per unit vol- 
ume, u> = fl/V, can thus be written as 



t 



(8) 



In the polydisperse limit, the p pt turn into a reservoir 
polymer density distribution p p (cr p ) and the sum over t 
into an integral over a p so that 



As- | da p p p (a p )e-^^ 



(9) 



We note that a similar expression was derived in [20] 
for the case of monodisperse colloids. This approach 
used fundamental measure theory and leads to the scaled 
particle-theory expressions for /hs and p^(a). Only 
colloidal gas-liquid demixing was investigated, and the 
polymer polydispersity was interpreted as arising from 
the compressibility of polymer chains of a given length. 
A semi-grand approach, with the slight modification of 
keeping the total polymer number fixed, is then appro- 
priate because the different polymer "species" can effec- 
tively be transformed into each other. 

For genuine polymer chain length polydispersity, how- 
ever, we are in practice dealing with a sample contain- 
ing a fixed number of polymers of each species. We 
therefore need to transform back to the canonical en- 
semble to obtain the Hclmholtz free energy density, / = 
u> + ^2 t p p tP P t- From (8) and bearing in mind the defini- 
tion of p pt , the polymer densities and chemical potentials 
are related by the following transformation: 



Ppt = - 



duj 



dp, 



V 



-Mc X ( CT pt) 



Ppt 1 



-fuf^Pt) 



(10) 



Simple algebra then gives for the free energy density 
/ = As + ^Pp t (lnp pt - 1) + ^p P *Mc X Kt) (11) 



Here we have separated off the ideal polymer part in 
the second term, dropping the contribution p pt In Ap t 
which is linear in the densities and therefore leaves the 
phase behaviour unaffected. In the limit of fully poly- 
disperse polymers the free energy (density) would then 
read 

/ = As + J do- p pp(crp) [lnpp(o-p) - 1] 

+ J da pPp (a p )p c h x s (a p ) (12) 

However, in the present study we neglect polydispersity 
of the polymers and focus on the effects of colloid poly- 
dispersity. After separating the free energy /h s of the 
pure colloid system into its ideal and excess parts, the 
total free energy (11) then simplifies to / = / ld + / ox , 
with 



/ ld = J da c p c (o- c )[\n p c (cr c ) - 1] +p p (lnp p 



-1) 



r 



rex i _ . ,ex/„. \ 
As +P P Mhs( cr p) 



(13) 



Importantly, evaluation of the free energy (13) requires 
as nontrivial input only the properties of the pure col- 
loid system, i.e. its excess free energy and the excess 
chemical potentials p^(cr) which are obtained by differ- 
entiation, Mhs( CT ) = Whs /$Pc( a )- To specify the free en- 
ergy fully, we therefore only need to assign appropriate 
expressions for fg£ in the colloidal fluid (or gas/liquid) 
and solid phases. For the fluid part of the excess free 
energy the most accurate approximation available is the 
generalisation by Salacuse and Stell [25] of the BMCSL 
equation of state [26, 27] while for the solid we adopt 
Bartlett's fit to simulation data for bidisperse hard sphere 
mixtures [28, 29]. Our previous work [1, 2] on polydis- 
perse hard spheres has shown that with these free en- 
ergy expressions quantitatively accurate fits to simula- 
tion data are obtained, and so we continue to use them 
for the present study. 

Some care is needed in the application of the excess 
free energy expression for the solid. The simulation data 
from which this is derived were obtained for mixtures 
of two species of hard spheres with diameters differing 
by no more than w 15%. The resulting excess chemi- 
cal potentials /x[! x (cr) are therefore accurate only within a 
small range around the mean colloid diameter. Outside 
this range they are unreliable; for example the limiting 
behaviour predicted from the Widom insertion principle 
for n — > is not retrieved correctly. This causes a diffi- 
culty because in the colloid-polymer interaction term in 
the free energy (13) we need the excess chemical potential 
evaluated at the polymer diameter, which we will take to 
be rather smaller than the mean colloid diameter. In fact, 
it is only in this regime of size ratios £ well below unity 
that the AO-model approximation of treating polymers 
as spherical particles is reasonable; for larger polymers 
the colloids can "see" the polymer chain structure [30]. 

To circumvent this problem, we follow previous work 
within the free volume approach [10, 17, 18] and always 
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evaluate the excess chemical potential p,° x (°~p) governing 
the polymer-colloid interaction from the BMCSL free en- 
ergy. The underlying physical approximation is that the 
free volume available to the polymers is not drastically 
different for colloidal fluid and solid phases so that the 
fluid (BMCSL) expression can also be used to estimate 
the free volume in colloidal solids. 



III. MOMENT FREE ENERGY METHOD 

A key feature of the excess free energies which we use 
to describe the colloidal fluid and solid phases is that they 
are truncatable [23]: they depend only on the finite set 
of moments Pi — J da c p c {a c )a % c (i = 0, . . . , 3) of the col- 
loidal density distribution p c (o~ c )- Here p c {a c )da c is the 
number density of colloids with diameters in the range 
<r c . . . <7 C + da c . The excess chemical potentials are then 
third order polynomials, 

The coefficients /i£ x , here are the moment excess chem- 
ical potentials of a pure hard sphere system, p,^ i — 
df^/dpi, and also depend only on the p,. The excess 
part of the free energy (13) of our interacting colloid- 
polymer mixture thus takes the form 

r x = /h c s x (W)+ppE^X- as) 

i 

Recall that in the last term of (15), i.e. in the polymer- 
colloid interaction, we always use the i derived from 
the BMCSL free energy. 

The excess free energy (15) depends only on a finite 
number of variables, namely the pi and the polymer 
density p p . Importantly, these can again be viewed as 
moments of an enlarged density distribution which col- 
lects all the densities of our system, namely (p c (<7 c ), p p ). 
Specifically, we define moments by 

Pi = J d<r c p c (a c )Wi(a c ) + ppwf (16) 

in terms of weights (w^(a c ) 7 wf) which are made up of 
a weight function w^(a c ) for the colloid part and a sin- 
gle coefficient wf for the polymer part. The colloidal 
moment densities, pi, (i = 0, . . . , 3) are then given by 
the weight functions (ct*,0) while the polymer density 
p p = pi is the moment with weight (0, 1). 

In summary, our system is described by an excess free 
energy which depends on five moment densities of the 
enlarged density distribution (p c (<7 c ), p p ) and is therefore 
by definition truncatable. This allows us to employ the 
moment free energy method for the calculation of phase 
diagrams. As described in [23, 31-33], and for the spe- 
cific case of polydisperse hard spheres in [2], this maps 
the full free energy (13), with its dependence on all details 



of (p c (cr c ), p p ) through the ideal part, onto a moment free 
energy (MFE) / mom ({pJ) = (£i V, - p tot ) + / ox (M) 
which depends only on the moments pi (i = 0, ... ,4). 
Here p to t — po + pi is the total number density of col- 
loids and polymers and the \ are Lagrange multipliers 
which depend implicitly on the values of the {pi}. From 
this MFE, phase behaviour can then be found by the 
conventional methods for finite mixtures, treating each 
of the pi as a number density of an appropriate quasi- 
species. For truncatable free energies this locates exactly 
the cloud points, i.e. the onset of phase separation, as 
well as the properties of the coexisting "shadow" phases 
that appear there. Inside the coexistence region, one in 
principle needs to solve a set of highly coupled nonlinear 
equations and the predictions derived from the MFE are 
only approximate. However, by retaining extra moments 
with adaptively chosen weight functions, increasingly ac- 
curate solutions can be obtained by iteration [23, 34]. 
Using these as initial points, we are then able to find full 
solutions of the exact - for our model free energy - phase 
equilibrium equations [35, 36]. 

We work with dimensionless units in the following. We 
call the overall colloid density distribution in the system 
the "parent" colloid distribution, pi°\a c ), with overall 
colloid density pc°' = J da c p^\a c ). The normalised 

parent size distribution is ni°\a c ) — pi°\cr c ) / p[ ^ ; we 
denote its mean diameter by cto- Lengths are measured 
in units of o~o; the dimensionless polymer diameter <r p 
then coincides with the polymer-to-(mean) colloid size 
ratio £. All densities are referred to the volume of a unit 
colloid particle, (7t/6)<7q. 

IV. PHASE DIAGRAM TOPOLOGIES 

In this section we will describe our results for the over- 
all phase behaviour of a mixture of polydisperse hard 
spheres and monodisperse polymers. Our numerical work 
requires a choice to be made for the colloidal parental 
diameter distribution. We focus on a triangular distribu- 
tion, 

(0) / \ 1 J cr — (1 — w) for 1 — w < <J C < 1 

n c (^J - ^ | (1 + w ) - for 1 < (t c < 1 + to 

whose width parameter w is related to the polydispcrsity 
by w = \/6 5. For the moderate values of 5 of interest 
here one expects other distribution shapes to give qual- 
itatively similar results. As in the case of polydisperse 
hard spheres without added polymer, this is based on the 
intuition that for narrow size distributions S is the key 
parameter controlling the phase behaviour [1, 2, 37]. 

To assess the range of possible phase diagram topolo- 
gies we consider initially four combinations resulting from 
having either small or large polymers, £ = 0.2 and 
0.4, and colloids with small or moderate polydispersity, 
S = 0.05 and 0.08. Fig. 2 shows the resulting phase di- 
agrams in the 4> c -4>p plane. The nature of the different 
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FIG. 2: Phase diagram on a 2 x 2 grid of values of (£, 6) € 
{0.2, 0.4} x {0.05, 0.08}. In each region the nature of the 
phase(s) coexisting at equilibrium is indicated (F: fluid, S: 
solid, G: gas, L: liquid). The dashed lines indicate the best 
guess of the phase boundary in regions where our numerical 
data become unreliable. 
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FIG. 3: Phase diagram for polydisperse hard spheres without 
added polymer. Note the appearance of regions with multiple 
solid coexistence as polydispersity increases. Horizontal cuts 
at 5 — 0.05 and 0.08 give the behaviour along the baselines of 
the graphs in the previous figure. From [1] . 



coexisting phases, i.e. gas, liquid, fluid and solid, is indi- 
cated in each region. Probably the most striking feature 
of these phase diagrams is the presence of multiple solids: 
in addition to the phases that are found for monodisperse 
colloids (see Fig. 1), size polydispersity in the colloids 
allows the system to phase separate into several solid 
phases. 

Qualitatively, the phase diagrams of Fig. 2 can be un- 
derstood by combining the behaviour observed in two 
simpler cases: the limit of a mixture of polymers and 
monodisperse colloids, which gives the phase diagrams 
sketched in Fig. 1, and the polymer- free limit of poly- 
disperse hard-sphere colloids [1, 2], with phase diagram 
shown in Fig. 3. The first limit clarifies why we need to 
distinguish between small and large polymer sizes: the 
short- or long-range nature of the depletion interaction 



dictates whether or not gas-liquid phase separation oc- 
curs. The second case of pure colloids clarifies the role of 
polydispersity in the phase behaviour. In fact the base- 
line (4> p = 0) of each diagram in Fig. 2 corresponds to a 
horizontal cut at the appropriate colloid polydispersity S 
in the polymer-free phase diagram of Fig. 3. The multi- 
ple solid phases occurring in Fig. 2 are thus "inherited" 
from the phase behaviour of the polydisperse colloid sys- 
tem without polymer. 

Having explained the broad intuition behind Fig. 2, let 
us look at the individual phase diagrams in more detail. 
Consider first the case with small polymers, £ = 0.2, and 
moderate colloid polydispersity, 5 = 0.05. From Fig. 3 
we see that as we move along a dilution line at 5 = 0.05, 
i.e. as we increase the colloid density at fixed normalised 
parent size distribution, we have the following sequence 
of phase splits: F -> F+S -» S -> S+S This se- 

quence re-occurs on the baseline of the bottom left graph 
of Fig. 2 as it must. If we now gradually add polymers 
to the solution we find that, as in the monodisperse case, 
the F+S coexistence region becomes wider; but in the 
polydisperse case, the boundary between the F+S and 
S regions eventually meets that between the S and S+S 
regions, resulting in a triple point that marks the be- 
ginning of an F+S+S coexistence region. A phase split 
of this type is not present in the polymer-free system 
at this polydispersity 5 and results from the interplay 
of the polymer-generated attraction force, which favours 
fluid-solid phase splits over single-phase solids, and the 
polydispersity of the colloids. At higher values of poly- 
mer or colloid concentration, coexistence of a fluid with 
an increasing number of solids then occurs. 

At higher colloid polydispersity, 5 = 0.08, the phase 



diagram topology is different. On the 



baseline 



the phase split sequence now consists of a fluid phase 
coexisting with an increasing number of solid phases, at 
least up to 4> c i=a 0.63. As Fig. 2 (top left) shows, all 
these phase boundaries are affected by the addition of 
polymer in the same way as the boundary between F 
and F+S in the monodisperse case, shifting to smaller 
4> c as <ft p is increased. There are no phase boundaries 
at which a fluid phase is lost and which would move to 
larger <ft c , as was the case for the F+S — > S boundary at 
S = 0.05. Consequently no phase boundaries meet with 
increasing <j> p and the phase split sequence remains as for 
the polymer-free system. 

Using similar arguments it is possible to explain the 
topologies at £ > £ c . In the system with polydispersity 
S = 0.05 shown in Fig. 2 (bottom right), we see again 
that the polymer-induced attraction favours fluid-solid 
splits. As the polymer concentration is increased from 
the baseline - which is the same as in Fig. 2 (bottom 
left) - the phase equilibria therefore acquire an additional 
fluid phase, via transitions from F to G+L, from F+S 
to G+L+S and from S to F+S. Because of the colloid 
polydispersity, however, the same mechanism now also 
causes coexistence of a fluid with several solids, with e.g. 
S+S becoming F+S+S. At higher level of polydispersities 
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FIG. 4: Phase diagram at £ = 0.3 and S = 0.08. Note the 
topology at high polymer concentrations where the region of 
gas- liquid phase splits terminates, giving way to phase separa- 
tion involving only a single fluid phase and one or more solids. 
The arrows indicate at which points in the phase diagram the 
diameter distributions in Fig. 13 below are calculated. 



(5 = 0.08, Fig. 2 (top right)), all phase splits on the 
baseline already contain a fluid phase which splits into 
G+L on adding polymer. In particular, this mechanism 
turns F+S+S into G+L+S+S, a polydispersity-induced 
coexistence of gas, liquid, and two solids. 

What remains unclear in this last phase diagram 
(Fig. 2 (top right)) is how the phase boundaries meet 
at high polymer concentration. In order to understand 
this we explored an intermediate polymer size £ = 0.3 
at the same polydispersity S = 0.08. This is easier com- 
putationally than exploring higher polymer densities at 
£ = 0.4 because the smaller value of £ — 0.3 moves the 
interesting region of the phase diagram to lower (f> p . As 
Fig. 4 demonstrates, the region involving gas-liquid sep- 
aration eventually terminates at high density, forming a 
closed "loop" ; above this loop one has a single fluid phase 
coexisting with one or multiple solid phases. 

The phase diagram topologies that we have shown so 
far are the most common that we have encountered, but 
are not the only ones possible. This is clear from two 
arguments. First, we have not exhausted all possible se- 
quences of phase splits along the polymer-free baseline; 
for colloid polydispersity <5 = 0.07, for example, Fig. 3 
shows that the sequence is F — > F+S — > F+S+S — > S+S 
— > S+S+S Second, at constant colloid size dis- 

tribution the phase diagrams need to change smoothly 
as the polymer size £ is varied, connecting e.g. the phase 
diagrams on the left of Fig. 2 to those on the right. In 
Fig. 5 we show a sketch of the various topologies that 
should result as colloid polydispersity S and polymer size 
£ are varied. The usual rule that exactly one phase must 
be lost or gained on crossing a boundary and the con- 
dition of smooth variation with £ and S constrain the 
possible topologies to a large extent. They do not, how- 
ever, make them fully unique, and Fig. 5 shows a pos- 



sible alternative topology for one of the phase diagrams 
(middle right). Nevertheless, we believe that the scenario 
shown in Fig. 5 is the most physically plausible sequence 
of phase diagram topologies. The topology shown in the 
top right corresponds to Fig. 4 and the top right of Fig. 2, 
although the numerically calculated phase diagrams do 
not extend sufficiently far to see all the different phase 
splits expected theoretically. The remaining phase dia- 
grams in Fig. 2 correspond similarly to top left and bot- 
tom left and right of Fig. 5. 

One intriguing feature of the predicted topologies is 
that phase splits involving gas and liquid phases can per- 
sist even after the region of pure G+L phase equilibrium 
has vanished from the phase diagram. This is visible in 
the middle and top rows of Fig. 5. As £ is decreased, 
the gas-liquid loop retracts relative to the onset of fluid- 
solid phase separation. At the value of £ where the two 
boundaries cross, the area of the G+L-region has shrunk 
to zero. However, a three-phase G+L+S-region of finite 
size survives at this point, and cannot disappear discon- 
tinuously as £ is varied. Thus phase equilibria involving 
gas and liquid phases can occur even if the initial phase 
separation that occurs on increasing density is always 
into fluid and solid phases. This effect is possible only 
because of the size polydispersity of the colloids. 



V. QUANTITATIVE ANALYSIS OF PHASE 
DIAGRAMS 

So far we have discussed the possible topologies of 
the phase diagrams of our colloid-polymer mixture. In 
the present section we will study their quantitative de- 
pendence on the colloid polydispersity as well as on the 
polymer-colloid size ratio £. We first analyse the be- 
haviour of the cloud and shadow curves which define the 
onset of phase coexistence, and consider the dependence 
of the crossover value £ c of the polymer size on colloid 
polydispersity. Fractionation effects are analysed next, 
both for fluid-solid and for gas-liquid phase separation. 
Finally we study the shape of the internal phase bound- 
aries between regions with two or more phases, showing 
how some general features can be predicted with simple 
arguments. 



A. Cloud and shadow curves 

The cloud curve defines the onset of phase separa- 
tion coming from a single-phase region, while the shadow 
curve records the properties of the incipient phase at this 
point. We now ask how these curves are affected quanti- 
tatively by colloid polydispersity. Looking back at Fig. 2, 
it is clear that at fixed £ any such effects on the cloud 
curves are small. In Fig. 6 we show the cloud curves 
together with the shadows at 5 = 0.05 and 5 = 0.08. 

Also the shadow curves are seen to be only weakly 
affected by polydispersity. The dotted lines connect- 
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FIG. 5: Sketch of expected phase diagram topologies as a function of polymer size £ and colloid polydispersity 6. The sequence 
of phases along the P = O-baseline in each row corresponds to 8 — 0.05, 0.07 and 0.08 respectively. Middle row, right: Dashed 
lines indicate a possible alternative topology which can still be connected smoothly to the one below but is physically less 
plausible. Bottom row, middle: The dotted line indicates the remnant of the G+L+S three-phase region, whose area shrinks 
to zero at the transition between the two topologies on the left and right. 



ing cloud and shadow pairs demonstrate that, as in the 
monodisperse case, there is strong colloid-polymer frac- 
tionation at the onset of phase separation, with the poly- 
mer preferably found in the lower-density phase. How- 
ever, polydispersity does introduce an additional feature, 
namely colloid size fractionation; this is discussed in the 
next subsection. We consider next the quantitative ef- 
fects of polydispersity and polymer size on the onset of 
gas- liquid phase separation as shown in Fig. 6 (right). 
To do so, we calculate the location of the gas-liquid crit- 
ical point and the associated gas-liquid spinodal. These 
quantities can be obtained in a straightforward way and 
without approximation, from the MFE. The conceptual 
reason for this is that the local stability of a system can 
be evaluated by only checking the effects of perturbations 
along the moments contained in the excess part of its free 
energy [23]. The results are shown in Fig. 7(a) where we 
plot the previously derived cloud curve at £ = 0.4 and 
8 = 0.05 together with spinodals and critical points ob- 
tained at different values of the polymer size £. The 
dependence on £ is pronounced. In Fig. 7(b) we contrast 
this with the change in the location of the critical point 
as 8 rather than £ is varied. The effect is very much 
smaller, and would be invisible on the scale of Fig. 7(a). 
We can nevertheless indicate the direction of the shift in 



the critical point with increasing 8. This is seen to point 
inwards, which means that colloid polydispersity tends 
to delay the onset of gas-liquid separation. 

We note briefly that the shift in the critical point coor- 
dinates 4> c and (f) p is quadratic in 5 for the modest poly- 
dispersities studied here. This arises because the critical 
point condition only involves moments of the size distri- 
bution, all of which are shifted by amounts of 0(5 2 ) for 
small 8. The dependence of the phase boundaries on £, on 
the other hand, is essentially identical to that observed 
for monodisperse colloids. The spinodals - and, by infer- 
ence, the cloud curves - shift to the right as £ decreases; 
eventually, at the threshold value £ — £ c , the G+L phase 
separation then becomes metastable with respect to the 
F+S transition. 

As reviewed in Sec. I, a number of studies have at- 
tempted to locate this threshold value £ c of the polymer 
size, leading to estimates of ~ 0.3... 0.6 in theoretical 
work and ss 0.25 in experiments. The effect of colloid 
polydispersity on £ c has remained unclear, however. To 
address this point, we have investigated different pairs 
of values for (£, 8) and determined the cloud curve for 
each of them. This is computationally nontrivial be- 
cause the equilibrium equations need to be solved very 
accurately to determine the point at which a fluid-solid 
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FIG. 6: Plots of cloud (thick line) and shadow (thin line) 
curves at 8 = 0.05 and 0.08. Dotted lines connect exam- 
ple cloud-shadow pairs. Left: small polymers, £ = 0.2; here 
the cloud curve marks the onset of fluid-solid phase separa- 
tion. Right: large polymers, £ = 0.4. The cloud curve at the 
larger <f> p gives the onset of gas-liquid separation. Because 
cloud and shadow curves at either side of the critical point 
(marked by the circle and triangle) are almost indistinguish- 
able, we plot only the cloud curve below the critical point and 
the corresponding shadow above it. At low <j>p °ne has onset 
of fluid-solid coexistence; the squares mark the triple points 
where the two branches of the cloud curve meet. 
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FIG. 7: (a) Spinodal curves obtained for different values of 
£ together with the cloud curve at (£ = 0.4, 8 = 0.05) and 
critical point, (b) Change in the position of the critical point 
with colloid polydispersity 8, at fixed polymer size £ = 0.4. 




binodal in a polydisperse system interchanges stability 
with a gas-liquid binodal. Fig. 8 summarises our results 
for two different colloid polydispersities. For 6 = 0.01, 
we see that the gas-liquid binodal becomes unstable at 
Cc = 0.31 ± 0.01, while at 5 = 0.05 this is shifted to the 
lower value £ c = 0.29 ± 0.01. We see that polydisper- 
sity favours gas-liquid over fluid-solid phase separation, 
shifting the threshold polymer size above which stable 
gas-liquid phase coexistence is observed to lower values. 
Indeed, for a system with larger colloid polydispersity 
(S = 0.1, data not shown) we find a further significant 
shift to£ c = 0.25±0.01. This is in good agreement with 
the experimental data [5] - though for a somewhat larger 
polydispersity than quoted in [5] - and is significantly be- 
low the value £ c = 0.32 found in [10] for monodisperse 
colloids. In summary, colloid polydispersity has signif- 
icant effects on the threshold value £ c for the polymer 
size, with a 5 of only 10% reducing it from £ c w 0.31 to 
£ c « 0.25. 

Several further comments are in order. First, for small 
5 one would expect £ c to be decreased from its monodis- 
perse value by terms of order <5 2 . This is because phase 
boundaries should generically shift by terms of this or- 
der [21, 22] as polydispersity is increased, while they 
should react smoothly, i.e. linearly, to changes in £. Un- 
fortunately, we do not have at present sufficient (and 
sufficiently precise) data to verify this expectation. Sec- 
ond, our predicted £ c = 0.31 for S = 0.01 is just slightly 
smaller than the monodisperse value £ c = 0.32 from [10]. 
This could be due to the small polydispersity in our 
case, though if shifts in £ c are indeed of order S 2 then 
S = 0.01 should give a shift 100 times smaller than that 
(0.31 - 0.25 = 0.06) observed for 5 = 0.1, which would 
be negligible. More likely the small difference is due to 
the fact that in [10] the scaled-particle or Percus-Yevick 
expression was used for the free volume term Mhs( a p) 
in (13) while we employ the BMCSL expression, consis- 
tently with the free energy fg£ describing the colloidal 
fluid. Finally, we recall from Sec. VA that colloid poly- 
dispersity delays the onset of gas-liquid coexistence. On 
the other hand, we just saw that it favours gas-liquid 
separation relative to fluid-solid. The conclusion is that 
polydispersity must disfavour fluid-solid coexistence, and 
do so more strongly than for gas-liquid. This is quite 
plausible, since polydispersity causes inefficient particle 
packing in a crystalline structure, while it can actually 
increase packing efficiency in a liquid. 



Fractionation effects 



FIG. 8: Plot of the phase boundaries defining the onset of 
F+S (solid line) and G+L (dashed line) coexistence, at colloid 
polydispersities 8 = 0.01 (left) and 8 — 0.05 (right). To make 
the curves at different polymer sizes £ visually distinguishable, 
those for the lowest and highest ^-values have been shifted 
along the horizontal c -axis by —0.1 and 0.1 respectively. 



As pointed out above, in colloid-polymer mixtures with 
polydisperse colloids the onset of phase separation causes 
fractionation not just between the colloids and the poly- 
mer, but also between the various sizes of colloid parti- 
cles. To demonstrate this, Fig. 9 (left) shows how the 
polydispersities of the fluid and solid phases vary along 
the cloud and shadow curves of Fig. 6 (left). On the x- 




FIG. 9: Left: Polydispersity 8 of the fluid cloud and solid 
shadow phases versus their colloid volume fraction <j> c , at the 
onset of F+S coexistence for £ = 0.2 and parent polydispersity 
0.05 and 0.08. Sec Fig. 6 (left) for the corresponding poly- 
mer volume fractions <fr p . Solid lines refer to the fluid cloud 
phase, dot-dashed lines to the solid shadow; dotted lines con- 
nect sample cloud-shadow pairs. Middle and right: Example 
plots of the normalised colloid size distributions of the phases 
indicated by the circles in the left plot. 



axis we plot the colloid volume fraction <j) c of the phases, 
so that the phases are now represented in terms of (cf> c , S) 
instead of (<j) c , <p p ) as before. The corresponding values of 
(j) p can of course be read off from the cloud and shadow 
curves in Fig. 6 (left). In particular, the ends of the 
lines in Fig. 9 (left), marked by empty circles, correspond 
to the limit </> p — > 0. Examples of the normalised col- 
loid size distributions, n c (er c ) = p c (a c )/ j da c p c (a c ), are 
shown in Fig. 9 (middle and right). These results give 
a concrete demonstration of the fact that the depletion 
interaction can be used to systematically control and re- 
duce the polydispersity of colloidal mixtures, as has been 
known for a number of years. In particular Bibette [38] 
suggested a procedure where colloid samples are fraction- 
ated from solutions of volume fraction <p c w 0.1. As we 
can see from Fig. 9 (left), this is indeed the regime where 
fractionation is strongest, i.e. where the shadow phase 
has a much smaller polydispersity than the parent (which 
is identical to the cloud phase). For phase separation at 
larger <j> c , corresponding from Fig. 6 (left) to smaller poly- 
mer volume fractions <fi p , the reduction in polydispersity 
is less pronounced. 

Colloidal size fractionation effects do of course occur 
not only in fluid-solid phase separation, but also in gas- 
liquid coexistence. Evans et al [21, 22] used a pertur- 
bative approach to study such fractionation effects for 
near-monodisperse parent size distributions. They pre- 
dicted that the difference in mean particle diameters of 
two coexisting phases, Aa c , should univer- 

sally scale as S v for small S, with exponent r] = 2. In [9] 
Fairhurst and Evans verified this relation experimentally 
for a colloid-polymer mixture: they used solutions of col- 
loidal PMMA with random polystyrene coils, with size 
ratio £ = 0.45. Collecting data from a number of sam- 
ples with parental colloid and polymer volume fractions 
in the region </> c = 0.15 . . . 0.5 and </> p = 0.15 . . . 0.3, they 
estimated a power-law exponent r\ = 2.16 ± 0.44. 

To compare our numerical results with these exper- 
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FIG. 10: Log-log plot of the difference in mean colloid diam- 
eter between coexisting gas and liquid phases, as a function 
of the parent polydispersity. The polymer- colloid size ratio 
is £ = 0.45. Solid and dashed lines show our predictions for 
two different choices of colloid and polymer concentration, as 
indicated in the legend. The dotted line is a power law with 
exponent ij = 2. The circles, together with grey error bars, 
indicate the experimental data points from [9], which were 
obtained from a set of samples in the range 4> c = 0.15 . . . 0.5, 
(p p = 0.15... 0.3. 



imental data, we considered the same polymer-colloid 
size ratio £ = 0.45 and two different choices for colloid 
and polymer concentrations, (p c = 0.3, <j) p = 0.2) and 
(p c = 0.2, 4> p = 0.3), with colloid polydispersity rang- 
ing from 6 = 0.025 to 0.40. We fix p c rather that <p c 
here since this is the case considered in the perturbative 
theory. The colloid volume fraction <f> c = (cr^)p c then 
increases with 6 and lies between 0.2 (for p c = 0.2 and 
small 5) and 0.45 (for p c = 0.3 and the largest 5). The 
overall range of variation of <f) c and (j> p is thus comparable 
to but somewhat smaller than in the experiments of [9]. 
Fig. 10 shows our numerical predictions for the difference 
in mean diameter of the coexisting gas and liquid as a 
function of parent polydispersity, together with the data 
from [9] . As expected our calculations are consistent with 
the universal quadratic scaling at small polydispersities. 
More importantly, they also show reasonable quantita- 
tive agreement with the experimentally measured values. 
More precise experiments covering a narrower range of 
colloid and polymer concentrations would obviously be 
useful, to permit a more stringent test of our calcula- 
tions. 



C. Inner phase boundaries 

Having discussed the properties of the cloud and 
shadow curves, and the fractionation effects which occur 
at or near the onset of phase separation, we now turn to 
the internal phase boundaries and their dependence on 
the polydispersity and polymer-colloid size ratio. 

A general trend apparent from Fig. 2 is that, with 
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increasing polymer concentration, the phase boundaries 
shift systematically to lower colloid volume fraction. 
Phase separation into several solids, for example, can 
thus occur at much lower colloid concentration than in 
a system without polymer. An intuitive explanation for 
this is that the osmotic pressure in the system is increased 
as polymers are added. This compresses the colloids and 
can be viewed as effectively increasing the colloid volume 
fraction. From Fig. 3, phase equilibria involving several 
solids will then occur earlier, as we observe. 

Quantitatively, perhaps the most striking feature of 
Fig. 2 (right) is the fact that for large polymers the 
internal phase boundaries which continue those in the 
polymer-free system are almost linear in the (</> c ,</> p )- 
plane. For the boundaries of the G+L+S region in Fig. 2 
(bottom right) this is straightforward to rationalise. For 
monodisperse colloids, the colloid-polymer mixture con- 
tains only two different particle species and so the bound- 
aries of the three-phase G+L+S triangle must be exactly 
straight; compare Fig. 1. Introducing a small degree of 
polydispersity S should then leave these boundaries ap- 
proximately straight as observed. This explanation al- 
ready falters for Fig. 2 (top right), however, where the 
G+L+S region widens rather than narrows as <fi p in- 
creases, precluding a straightforward analogy with the 
three-phase triangle in a monodisperse system. It fails 
completely for the inner phase boundaries marking the 
transition between phases involving several solids, which 
are induced purely by polydispersity effects and have no 
monodisperse analogue. A different explanation is there- 
fore required: we will see that the near-linearity of the 
inner phase boundaries arises from the fact that the gas 
phase consists almost entirely of polymers, plus a negligi- 
ble amount of colloids. It can therefore act as a "polymer 
buffer" for the other phases, which contain only small 
amounts of polymer and are otherwise similar to those in 
the pure colloid system. For concreteness we focus below 
on a case where these colloid-rich phases are solids, but 
the argument applies equally well when there is a liquid 
among them. 

To substantiate this hypothesis, we have extracted in 
Fig. 11 the properties of the coexisting phases for poly- 
mer size £ = 0.4 and polydispersity S — 0.05, along the 
phase boundary between the F+S+S and F+S+S+S re- 
gions highlighted by the bold line in Fig. 11 (a). The 
position along the boundary in the plots (b-f) is pa- 

rameterised by the parent polymer volume fraction p ^ , 
i.e. the y-coordinate from (a). In (d) and (c) we show 
the polymer and colloid volume fractions in the various 
phases. We see that, as anticipated, the fluid phase con- 
sists almost entirely of polymers, i.e. it is an extremely 
dilute colloidal gas. The solids, on the other hand, are es- 
sentially pure colloid phases and contain almost no poly- 
mer. In addition, as plots (e) and (f) show, their proper- 
ties remain unchanged along the phase boundary and can 
be inferred from the extrapolation to the polymer-free 
limit </>p°' = 0. This applies in particular to their com- 
mon pressure IT* , and the average colloid volume fraction 
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FIG. 11: Properties of the coexisting phases, for a system 
with £ = 0.4 and S — 0.05, along the phase boundary shown 
in (a) by the dark line, (b) Fractions of system volume 
occupied by the various phases; the newly forming solid has 
vanishing fractional volume, (c) Fractional volumes of the 
solids normalised by the total fractional volume occupied by 
solid phases. (d,e,f) Polymer and colloid volume fractions and 
colloid polydispersity of the coexisting phases. 




0.58 O.f, 0.62 0.64 0.66 0.68 



FIG. 12: Left: Pressure plot of a polymer-free system with 
polydispersity S = 0.05. The values of the pressure at the 
phase transitions are marked by the horizontal lines. Right: 
Phase diagram of a colloid-polymer mixture with the same 
colloid polydispersity and polymer size £ = 0.4, plotted as p p ^ 
vs (j)^ . Two phase boundaries are extrapolated (dashed) to 
4>c^ = 0. The agreement between the extrapolated values of 
the intercepts, 13.7 and 18.7, and the equilibrium coexistence 
pressures II* in the polymer-free system, 13.0 and 18.4, is 
good. 



in the solids, <j>*. The fractional system volumes oc- 
cupied by the solids also need to remain in constant pro- 
portion to each other, to maintain the overall colloid size 
distribution, consistent with the results shown in Fig. 11 
(b). Figure 11 (c) shows, more explicitly, that when we 
normalise the fractional volumes of the solid phases by 
the total fractional volume of all solids, they become con- 
stant along the phase boundary. The normalisation fac- 
tor is J2 a ^ g v^ = 1 — v( g \ where "g" denotes the gas 
phase. 

We can now infer the shape of the phase boundary. 



Since only the solid phases contain significant concen- 
trations of colloid, particle number conservation requires 
that the parental colloid volume fraction is given by 

^(0) = (i (17) 

Similarly, because the polymers are found almost exclu- 
sively in the gas, we have 

«£(°> = w («)^(g) = v (g){3 p W ( 18) 

Here we have used the fact that in our density units the 
polymer density p p and volume fraction p are related 
simply by a factor £ 3 . In the same units, the pressure 
of the gas phase is Il( g ) = p p s ^ because the polymers are 
ideal and the colloids make only a negligible contribu- 
tion. At equilibrium this pressure must equal that in the 
colloid phases, II*, hence 

pM = IT (19) 

Combining the last three equations now gives the desired 
relation between the parent polymer and colloid volume 
fractions along the phase boundary: 

p (0) = r 3 (O) = n , ^ _ 4jjpj {2Q) 

This is indeed linear, as we set out to show, and obeys 
the general trend that with increasing <f) p ^ phase transi- 
tions occur at smaller parental colloid volume fractions 
4>i°\ The relation (20) also predicts that if we extrap- 
olate the straight phase boundaries to 0c°' = 0, i.e. to 
their intersection with the y-axis, the intercept will be 
4 0) = £ 3 IT, or p p 0) = IT if we plot p p 0) on the y-axis. 
In Fig. 12 we check this prediction explicitly for two of the 
numerically calculated phase boundaries of the system in 
Fig. 11, finding good agreement. More generally, one con- 
cludes from (20) that all the straight phase boundaries in 
systems with sufficiently large polymers should become 
approximately independent of the actual polymer size £, 
once they are plotted in terms of p p ^ versus 0c '' . The re- 
striction here is that the polymers must not be too close 
to the threshold size £ c - This ensures that the straight 
phase boundaries extend close to the polymer-free base- 
line, thus justifying our extrapolation to <p p = to deter- 
mine IT* and <fi*. For £ close to £ c , on the other hand, as 
in e.g. Fig. 4 above, the relevant polymer-free "reference 
system" can no longer be obtained by straightforward 
linear extrapolation. Nevertheless, our arguments show 
that along straight internal phase boundaries, and indeed 
along appropriate straight lines between the phase bound- 
aries, the colloidal size distribution of the dense phases 
should remain approximately unchanged. 

To illustrate this phenomenon, we plot in Fig. 13 three 
examples of normalised colloid diameter distributions at 
the three points marked in the phase diagram of Fig. 4. 
These points lie approximately on a straight line fol- 
lowing the direction of the phase boundaries. As ex- 
pected, we find that the size distributions in the liquid 
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FIG. 13: Normalised colloidal size distributions n c (a c ) in co- 
existing phases; the parent distribution is shown for compar- 
ison. The three graphs correspond to the values of (<j> c , (f> p ) 
indicated in Fig. 4. 



and solid phases remain approximately constant. This 
seems rather remarkable, given that the overall colloid 
and polymer fractions in the three cases are very differ- 
ent. 



VI. CONCLUSION AND OUTLOOK 

We have studied theoretically, for the first time, the 
equilibrium phase behaviour of mixtures of polydisperse 
hard sphere colloids and monodisperse polymers. Our 
treatment is based on the AO model, which treats the 
polymers as spherical particles interacting only with the 
colloids. Within the van der Waals approximation of 
Lekkcrkcrkcr et al, the polymer-colloid interaction is de- 
termined by the excess chemical potentials of the polydis- 
perse hard sphere system. We have therefore used as in- 
put suitable free energy approximations for the colloidal 
fluid and solid, choosing the BMCSL and Bartlett's "ge- 
ometric" free energy, respectively; these were shown in 
previous work to give quantitatively reliable predictions. 
Complete phase diagrams, taking full account of frac- 
tionation effects, are found exactly - to within numerical 
accuracy - by the use of the moment free energy with 
extra adaptive moments; this is possible because the free 
energies involved are truncatable. The intricate features 
of phase separation in pure polydisperse colloids, includ- 
ing fractionation into several solids, combine with the 
appearance of polymer-induced gas-liquid coexistence to 
give a rich variety of phase diagram topologies. 

We studied in some detail the dependence of these 
phase diagrams on colloid size polydispersity S and the 
polymer-colloid size ratio £. Even for the moderate values 
of S we consider (S < 0.10), polydispersity has a signifi- 
cant influence on the number of phases and the topology 
of the phase boundaries. This influence arises because the 
sequence of phase transitions that is observed in a "base- 
line" system of pure colloids changes significantly with 5. 
Starting from this baseline the addition of polymer then 
causes fluid phases to phase separate into colloidal gas 
and liquid, or a new fluid phase to appear if only solids 
were present in the polymer-free baseline system. The 
strength of this effect grows with polymer size £; this is 
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reasonable because both the range and strength of the ef- 
fective colloid-colloid interaction caused by the presence 
of polymer grows with £. The strength of this attraction 
also grows with polymer concentration, and consequently 
all phase boundaries shift to lower colloid concentrations 
when more polymer is added. In particular, we predict 
that phase splits involving two or more fractionated solids 
should occur at fairly moderate colloid volume fractions, 
where the phase separation kinetics should be substan- 
tially faster than for the corresponding higher volume 
fractions in a system without added polymer. This sug- 
gests the exciting possibility that such solid-solid phase 
splits, which have not hitherto been seen in experiments, 
may be observable in this regime. We note, however, 
that the growth kinetics of polydisperse crystals [39] may 
still cause deviations from the predicted equilibrium be- 
haviour even if the regime of multiple solids is kinetically 
accessible. 

We deduced the possible phase diagram topologies 
from the requirement that they need to change smoothly 
with both 6 and £. One intriguing prediction of this 
is that, in the presence of colloid polydispersity, stable 
three-phase gas-liquid-solid coexistence can persist with- 
out a corresponding two-phase gas-liquid region in the 
phase diagram; see the middle and top rows of Fig. 5. 
This effect should also be observable experimentally, for 
values of £ just below the threshold £ c at which the stable 
gas-liquid coexistence disappears. 

In our quantitative analysis of the phase behaviour, we 
found that this threshold value £ c depends significantly 
on colloid polydispersity, decreasing from £ c w 0.31 for 
S = 0.01 to £ c w 0.25 for S = 0.1. We saw that poly- 
dispersity delays the onset of both gas-liquid and fluid- 
solid separation, but that the effect is stronger for the 
latter, causing the observed shifts in £ c . We further con- 
sidered fractionation effects and found that in fluid-solid 
coexistence these are most marked for high polymer and 
low colloid concentration, in agreement with experimen- 



tal protocols that are used to reduce polydispersity in 
colloidal systems. Also for gas-liquid coexistence the cal- 
culated fractionation effects are in reasonable agreement 
with the experimental data, and follow the universal frac- 
tionation relation in the limit of small polydispersity. Fi- 
nally, we observed that, for sufficiently large polymer 
sizes £, boundaries between regions of the phase dia- 
gram containing two or more phases are close to straight. 
This is due to the occurrence of a dilute colloidal gas 
phase which acts effectively as a polymer buffer. We 
derived from this an estimate of the location and slope 
of the phase boundaries in terms of the properties of a 
corresponding polydisperse colloid system without added 
polymer. 

An interesting extension of this research would be to 
the case where the polymers are polydisperse - with a 
fixed parent density distribution - while the colloids are 
monodisperse, or even further to the scenario where both 
colloids and polymers are polydisperse. Work in this 
direction is in progress. One would also like to assess 
the effects of polymer non-ideality, for which a num- 
ber of theoretical approaches have been developed, based 
on perturbation theory around the 0-point [40], integral 
equations [41] and effective colloid-colloid interactions de- 
rived from simulations of self-avoiding walk (SAW) poly- 
mers [42] . The method that fits most naturally into our 
framework is to add second-order virial terms to account 
for polymer-polymer interactions, as done in [40]. How- 
ever, recent simulation studies [43] indicate that the AO 
model gives quantitatively reasonable results also for in- 
teracting (monodisperse SAW) polymers up to quite large 
sizes £ < 0.34. For larger £, where significant deviations 
arise [43], one has to account both for polymer interac- 
tions and also increasingly for the detailed polymer chain 
structure. An "AO + second virial" model which only 
addresses the first effect is therefore unlikely to be useful 
over a significant range of £. 
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